Orbital magnetic moment and extrinsic spin Hall effect for iron impurity in gold 
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We report electronic structure calculations of an iron impurity in gold host. The spin, orbital and 
dipole magnetic moments were investigated using the LDA+(7 correlated band theory. We show 
that the around-mean-field-LDA+U reproduces the XMCD experimental data well and does not 
lead to formation of a large orbital moment on the Fe atom. Furthermore, exact diagonalization of 
the multi-orbital Anderson impurity model with the full Coulomb interaction matrix and the spin- 
orbit coupling is performed in order to estimate the spin Hall angle. The obtained value 7s ~ 0.025 
suggests that there is no giant extrinsic spin Hall effect due to scattering on iron impurities in gold. 



PACS numbers: 72.25.Ba,71.70.Ej,71.15.Rf,85.75.-d 

During the last several years, a broad interest and at- 
tention have been devoted to the spin Hall effect (SHE) in 
semiconductors [1 and metals 2j. This effect amounts 
to an observation of a transversal spin current when a 
charge current is flowing through a solid. The SHE is 
caused by the spin-orbit coupling (SOC) and can occur 
even in non-magnetic solids [2|. 

Recently, an experimental observation of a giant SHE 
in Au/FePt has been reported _3i. The spin Hall con- 
ductivity of ^ 10^ n^^cm^^ and the spin Hall angle as 
large as ~ 0.1 were measured [3]. Guo et al. suggested 
the effect to be of extrinsic origin due to the Fe and Ft 
impurities in gold [4 . They reported local (spin) den- 
sity approximation (LSDA) plus Coulomb U (LDA-|-[/) 
solution for Fe in Au with a very large orbital magnetic 
moment ^ l.BfiB- 

The results of Ref. contradict the LDA+U calcu- 
lations presented in Ref. [F that reported a tiny Ml ~ 
0.02/iB. The value of Ml from Ref. [4] is clearly in- 
consistent with the experimental x-ray magnetic circular 
dichroism (XMCD) data for the ratio of M^ and the ef- 
fective spin moment Ms, Rls = 0.034 H]. Assuming 
Ms ~ 3/i_B leads to Ml ~ O.l^s which is of an order 
of magnitude smaller than the prediction of Ref. |4j . It 
was suggested [7] that the discrepancy between Refs. [4] 
and |S] are due to different choices of the Coulomb U. A 
large value of Ml is calculated with [/ = 5 eV [Ij , while 
much smaller ones are obtained with U — 3 eV [S] . 

In this work we revisit the electronic and magnetic 
structure of the Fe impurity in Au. We examine different 
flavors of the rotationally-invariant LDA+U method [8]: 
the "fully localized limit" (FLL) as well as the "around 
mean field" (AMF) version. The results for the orbital 
magnetic moment Ml are compared with the available 
experimental data. 

Both LSDA and LBA+U methods yield broken- 
symmetry static mean-field solutions with ordered spin 
and orbital moments, whereas the true dynamical so- 
lution of an impurity in a non-magnetic host exhibits 



Ms = 2{Sz) — and Ml ~ {L^) = when no external 
magnetic field is applied and no preferential direction for 
the orientation of the moments exists. In order to go 
beyond the static mean-field and to incorporate the dy- 
namical electron correlations, we employ the exact diag- 
onalization (ED) method to solve a multi-orbital single 
impurity Anderson model (SIAM) [9] whose parameters 
are extracted from LDA calculations. We evaluate the 
spectral density at the Fe impurity in Au and estimate 
the spin Hall angle due to skew scattering on the impu- 
rity. A relation between the electronic structure and the 
extrinsic SHE is discussed. 

As a computational model we use a FeAuis supercell 
chosen to keep Fe and its 12 nearest Au neighbors sepa- 
rated from other impurity atoms. No relaxation is per- 
formed as it is not essential for the closed packed fee 
structure. We use the lattice constant of elemental Au, 
a = 7.71 a.u. All calculations are performed making 
use of a relativistic version (with SOC) of LDA+U im- 
plemented in the linearized augmented plane wave (FP- 
LAPW) basis [101 . The radii of the atomic muffin-tin 
(MT) spheres are set to 2.3 a.u. (Fe) and 2.5 a.u. (An). 
The parameter R x i^max — 7.6 determined the basis 
set size and the Brillouin zone was sampled with 343 
k points. We checked that a finer sampling with 729 
k points does not modify the results. 

First, we apply the conventional LSDA with the von 
Barth and Hedin [11] exchange-correlation potential im- 
plemented within the relativistic FP-LAPW method [H] ■ 
Our results for the spin and orbital moments inside the Fe 
MT sphere, Ms and Ml, are compared with the results 
of other calculations in Table |l] In spite of a relatively 
small (16 atoms only) and unrelaxed supercell, present 
results are in a fair agreement with VASP [5 results for 
a substantially larger and relaxed supercell containing 
108 atoms, with the tight-binding LMTO (TB-LMTO) 
results for 55 atom supercell [13], and with the KKR- 
ASA calculations [14]. All calculations indicate a small 
value of Ml for Fe impurity in Au, which is typical for 3d 
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transitional metals and alloys. The calculated Ml /Ms 
ratio is substantially smaller than the experimental value 
Rls — 0.034 [6]. Typically, Ml is underestimated in 
LSDA due to the lack of orbital polarization. This leads 
to a smaller ratio Ml /Ms and explains the disagreement 
with experiment. 



TABLE 1. LSDA magnetic moments on Fe in Au (in units of 
Ms)- 
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FIG. 1. (Color online) The spin-resolved d-orbital DOS for 
Fe impurity in Au calculated with LSDA (top) and AMF- 
LDA+f/, (7 = 3 eV (bottom 
projected DOS for Fe. 
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The calculated c?-orbital density of states (DOS) for 
Fe atom and the first nearest neighbor Au atoms are 
shown in Fig.[l]^top). LSDA yields fully occupied Fe spin- 
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t states that are hybridized with shallow Au d states. 
The Fe spin-J, states are substantially more localized. 
Analysis of the projected DOS shows that the eg-like 
{dx2-y2 + d^z'^_r2) states become practically fully spin- 
polarized while the t2g-like {dxy + dzx 
only partially polarized, see Fig. [ijtop) 

Now we turn to the LDA-|-[/ calculations. We com- 
pare FLL and AMF variants of the rotationally-invariant 
LDA-fC/ method. The full local occupation matrix with 
all spin off-diagonal components is preserved. The dou- 
ble counting of the non-spherical d-states contributions to 
the LSDA and the LDA-|-t/ parts of the potential is cor- 
rected. The exchange J — 0.9 eV was used for Fe (Slater 
integrals F2 = 7.75 eV, = 4.85 eV). The Coulomb U 
was varied from 3 eV to 5 eV. 

The spin Ms, orbital Ml and dipole Md [inj 3d mag- 
netic moments are given in Table |ll] together with the 
occupation of the Fe atom d orbitals, Ud- Both FLL and 
AMF flavors of LDA+C/ lead to an enhancement of Ml 
with respect to the LSDA estimate. It is due to non- 
spherical Coulomb and exchange interactions which are 
incorporated in LDA-I-C/ [TO] and cause an additional or- 
bital polarization to that induced by the spin-orbit cou- 
pling. The value of Ml increases with increase of the 
Coulomb U . It is observed that the FLL double count- 
ing yields a substantially stronger enhancement of AIl 
than the AMF method. 

There is also a substantial magnetic dipole moment 
Md formed on the Fe impurity. When the spin-orbit 
coupling is included and spin polarization is allowed, the 
initial cubic symmetry is broken and only the tetragonal 
symmetry remains. This effect is rather small in LSDA. 
It becomes substantially enhanced in LDA+[/ due to the 
additional orbital polarization. This effect is visible on 
the AMF-LDA-Ft/ DOS shown in Fig. [ijbottom). The 
main difference between LSDA and LDA+C/ occurs in the 
spin- J, channel for the t25-like states; the dxy state peels 
off from the dzx and dzy states and becomes occupied. 

Experimental XMCD data are available for Fe impu- 
rity in Au m. The measured value for Rls = Ml/[Ms + 
7 Mo] = 0.034 is in a very good agreement with our AMF- 
LDA-f [/ results for U in the range between 3 eV and 4 
eV. On the basis of these calculations we conclude that 
a reasonable value of the Coulomb U for Fe impurity in 
Au host is « 3 eV. 

Our FLL results for U = 5 eV are fairly close to those 
of Ref. |4j where the LDA+U double counting was not 
specified. In this case, the calculated Rls — 0.21 ex- 
ceeds the experimental XMCD value by almost an order 
of magnitude. Therefore, the FLL-LDA-|-[/ method does 
not satisfactorily describe the electronic structure of Fe 
impurity in Au. 

Both the LSDA and LDA-I-C/ methods yield broken- 
symmetry mean-field solutions with non-zero M5 and 
Ml ■ This is because the part of the Coulomb interaction 
treated in the Hartree-Fock-like approximation is trans- 
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TABLE II. Magnetic moments (in ^b) and 3d occupation na 
of the Fe impurity in An host as a function of Coulomb U 



FeAuis 




FLL 




AMF 


U (eV) 


3 


4 


5 


3 4 


Ms 


3.18 


3.21 


3.29 


2.94 2.90 


Ml 


1.24 


1.36 


1.44 


0.16 0.22 


7Md 


2.36 


2.71 


3.57 


2.16 2.35 


Rls 


0.23 


0.23 


0.21 


0.03 0.04 


rid 


6.00 


6.00 


5.97 


6.03 6.03 



formed into the exchange splitting field. This exchange 
field is of the order of a few eV (see Fig. [T]) and by far 
exceeds any imaginable external magnetic field. Thus, 
the LDA+C/ method, most probably, provides a reason- 
able description of the local-moment systems in (strong) 
external magnetic fields. 

When no external magnetic field is applied and no pref- 
erential direction for the orientation of the moments ex- 
ists, neither LSDA nor LDA+C/ suffice. Recently, an 
attempt has been made to go beyond the static mean- 
field approximation and to solve the SIAM for the Fe im- 
purity in Au employing the Hirsch-Fye quantum Monte 
Carlo method [7]. The authors used a simplified three- 
orbital model with a diagonal Coulomb vertex and a spin- 
diagonal spin-orbit coupling only. These simplifications 
make an estimate of the accuracy of the quantitative re- 
sults reported in Ref. f7| difficult. 

In order to deal with the electronic structure of the Fe 
impurity in the absence of the external magnetic field, 
we apply the finite-temperature ED method [T71 to the 
complete five-orbital d shell subject to the full spherically 
symmetric Coulomb interaction, spin-orbit coupling and 
a cubic crystal field. The effective multi-orbital impurity 
Hamiltonian can be written as 

kma mcr 
mm' aa' 

+ E iykd^mah^a + h.C.) (1) 
kniG 

H" 2 ^ ^ Ufnm'm" m'" ^ma^m' a' ^m'" a' ^m" ai 

inm' m" 
m'" aa' 

where creates an electron in the d shell and fol. 

creates an electron in the "bath" which models those 
host-band states that hybridize with the impurity d shell. 
The bath is predominantly composed of s and p bands of 
Au. The impurity-level position ea and the bath energies 
Efc are measured from the chemical potential ^. Param- 
eters ^ and AcF specify the strength of the spin-orbit 
coupling and the size of the cubic crystal field on the im- 



purity. They are determined from LDA calculations as 
^ = 60 meV and Acf = 32 meV. The hybridization pa- 
rameters Vfe and the bath energies do not depend on 
m and a, and thus the Hamiltonian preserves orbital and 
spin angular momenta. This is a good approximation for 
Fe in Au, since the lower-symmetry components of the 
hybridization turn out to be considerably smaller than f 
and Acf- 

For the ED method to be applicable, the continuum 
of the bath states is discretized. The parameters e^, et 
and Vfc are chosen so that the impurity Green's function 
corresponding to the discretized Eq. ([T]) with U — Q ap- 
proximates the impurity Green's function from the LDA 
as closely as possible. Namely, we require several lowest 
moments of the respective densities of states to coincide, 
m(™) = m(^°^\ where M„ = / e»go^(e)de/ / go^(e)de 
|18) . The integrals run over a 1 eV wide interval cen- 
tered at the Fermi level, which confines the LDA impu- 
rity resonance. The actual values of the bath parame- 
ters are e^^ = 80 meV and vj^^^ = 220 meV when the 
index k is restricted to a single value and the bath con- 
tains 10 spinorbitals (bath L "d-l-lO spinorbitals"). For a 
bath twice as large we get e^^^' G {—310, 340} meV and 

VjI"^ e {140, 170} meV (bath IL "d-f-20 spinorbitals"). 
The position of the impurity level obtained from this 
procedure is subsequently shifted by a Hartree-like con- 
tribution in order to maintain the LDA impurity occu- 
pation Ud = 6.18 when the local Coulomb term is intro- 
duced. 
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FIG. 2. (Color online) d-electron spectral function of the im- 
purity model of Eq. ([TJ with [/ = 3 eV and two variants of 
discrete bath: 10 bath spinorbitals (top) and 20 bath spinor- 
bitals (bottom). 



After the parameters of the discrete impurity model 
are set, the band Lanczos method [IHl [2D] is utilized 
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to determine the lowest lying eigenstates of the many- 
body Hamiltonian and to calculate one-particle Green's 
function Ggj^jyi- The resulting d-orbital spectral func- 
tion Im(Ggj^]y[)/7r is shown in Fig. [2] for the two models 
of the bath and for the Coulomb interaction parameters 
[/ = Fo = 3 cV and J = 0.9 cV (F2 = 7.75 eV, F4 = 4.85 
eV). The inverse temperature /3 = 500 eV^^ was used in 
these calculations. Although the details of the spectral 
peaks depend somewhat on the particular choice of the 
bath, the overall structure of the spectrum with peak(s) 
in the vicinity of the Fermi level is preserved when the 
bath parameters are varied. The spin S — 1.91, orbital 
L = 2.21 and total J = 3.87 moments are calculated for 
the d shell from the expectation values {X'^) = X{X + 1), 
X = S, L, J. Individual components of the moments, 
(Sz) and (Lz), vanish so that the spin-orbital symmetry 
is preserved and neither spin nor orbital polarization is 
induced in the absence of the external magnetic field. 

Now we estimate the spin Hall angle from the skew 
scattering on the impurity with a local magnetic moment. 
Following Refs. [11111], we evaluate the spin Hall angle as 



7s = 



12(5i (cos 2(5n — cos 260 



25 



15 cos 2(5^ 



10 cos 2dn 



(2) 



where Si is the p-wave phase shift which is assumed 
to be small for the non-resonant scattering, \Si\ = 0.1. 
The d-wave phase shifts 62 for j = 5/2 and 62 for 
j — S/2 are related to the occupations and of 
the corresponding 3d subshells via the Friedel sum rule 



j/{2j -I- 1) [H]- The Hall angle 75 vanishes when 



all d orbitals are equally occupied {62 — S2), it grows 
as an increasing spin-orbit coupling favors the occupa- 
tion of the j — 3/2 subshell, and it eventually reaches a 
maximum -•y^™^''-* = 4(5i/5 when the j = 3/2 subshell is 
completely filled (71.3/2 = 4 and 62 = n). 

The occupation numbers Uj and the Hall angle 75 ob- 
tained for the Hamiltonian of Eq. ([!]) are listed in Ta- 
ble Hill for the two bath models introduced earlier and for 
U = 3 eV and 5 eV. Results of the non-magnetic LDA 
calculation and of an atomic-like calculation without any 
bath orbitals are shown for comparison. The angle 75 
increases when the Coulomb U is added, and keeps grow- 
ing with further increase of U. For a fixed value of U, 
the Hall angle decreases with increasing hybridization V 
since the spin-orbit splitting in the host band is negli- 
gible and the hybridization thus effectively reduces the 
spin-orbit effects in the Fe d shell. 

Our results are consistent with the measurements of 
Fert et al. [32] of the anomalous Hall coefficient ^ 0.01 
in dilute 3d noble metal alloys. The angle 75 « 0.025 we 
obtain is 50% smaller than the earlier theoretical estimate 
7s = 0.055 by Gu et al. ^ and substantially smaller than 
the "giant" 75 = 0.11 reported by Seki et al. [3 . Note 
that the giant SHE interpretation of the experimental 
results in Ref. [3j has been recently challenged also from 



the experimental viewpoint 



TABLE III. Impurity occupations Uj and the spin Hall angle 
7s obtained with two different bath models for U — 3 eV 
and 5 eV. Non-magnetic LDA calculation and an atomic-like 
calculation are shown for comparison. 



Model U (eV) Ud n^/j »-5/2 7s 

LDA - 6.18 2.62 3.55 0.008 

bath I 3 6.18 2.68 3.50 0.011 

bath I 5 6.18 2.85 3.33 0.021 

bath II 3 6.18 2.94 3.24 0.026 

bath II 5 6.18 2.94 3.24 0.026 

no bath 3-5 6.18 2.98 3.19 0.029 



To summarize, our calculations show that the AMF- 
LDA+U method with the Coulomb U around 3 eV re- 
produces very well the XMCD experimental data for Fe 
impurity in An host. The calculated orbital moment at 
the Fe atom, Ml = 0.16/iB, is almost ten times smaller 
than that reported by Guo et al. U . We explicitly show 
that the reason for this difference is not only in the use 
of a smaller value of J7 ffj, but also in the appropriate 
choice of the LDA+t/ flavor. Furthermore, using the ex- 
act diagonalization of a multi-orbital impurity model, we 
estimate the spin Hall angle due to the scattering on the 
Fe impurity in the Au host as 75 ~ 0.025. It is substan- 
tially smaller than 75 — 0.11 reported by Seki et al. [3]. 
We conclude that scattering off Fe impurities in Au does 
not yield a giant SHE. 
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